Polymers in crowded environment under stretching force: globule-coil transitions 
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We study flexible polymer macro-molecules in a crowded (porous) environment, modelling them 
as self-attracting self-avoiding walks on site-diluted percolative lattices in space dimensions d = 2, 3. 
The influence of stretching force on the polymer folding and properties of globule-coil transitions are 
analyzed. Applying the pruned-enriched Rosenbluth chain-growth method (PERM), we estimate 
the transition temperature Te between collapsed and extended polymer configurations and construct 
the phase diagrams of the globule-coil coexistence when varying temperature and stretching force. 
The transition to a completely stretched state, caused by applying force, is discussed as well. 
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I. INTRODUCTION 

Long flexible polymer macromolecules in a good sol- 
vent possess configurational statistics, which is perfectly 
captured by the model of self-avoiding random walks 
(SAW) on a regular lattice This corresponds to the 
regime of polymer coils which holds provided that the 
temperature T is above the so-called O-temperature. In 
this regime, the mean end-to-end distance of an A^-step 
chain scales as Rn ~ N" SAVJ , where i^saw > 1/2 is an 
universal exponent, which depends on space dimension 
d only. With lowering the temperature, the effect of 
monomer-monomer attraction grows and the polymer ra- 
dius shrinks. At T — Tq the effective repulsion due to 
the volume exclusion constraint is exactly balanced by 
attractive interactions. At this particular temperature 
a crossover occurs from high temperature SAW behav- 
ior to O-statistics with exponent v®(d = 2) = 4/7 |2|, 
v&{d > 3) = 1/2. Below the 9-temperature, the entropic 
effects, which make the polymer chain swell, are over- 
come by interaction energy and a collapse to the globule 
regime (with collapsed size exponent v c = 1/d) occurs. It 
is generally accepted [l[ that the coil-globule transition 
is of second order, in the sense that the density of an 
infinite globule is zero at T = Tq and increases contin- 
uously when T is lowered further; more presicely, it is a 
tricritical point with the upper critical dimension d c = 3. 

The coil-globule transition is one of the fundamental 
problems in polymer physics [H, [3, 0| , being deeply con- 
nected with problems like protein folding, DNA conden- 
sation and chromatin organization. The properties of 
polymers in the vicinity of the 0-point can be success- 
fully studied on the basis of self-attractive self-avoiding 
walks (SASAW), where a nearest-neighbor interaction is 
included. Estimates of the transition temperature Tq 
for flexible polymers have been so far the subject of nu- 



merous studies 0, H, 0, i, i, E3, [LI El. Note that the 
0-temperature is a non-universal quantity, depending on 
the lattice type; in particular it appears to be propor- 
tional to the connectivity constant fi of a SAW on a given 
type of lattice (recent estimates of n and 0-temperatures 
on simple square and cubic lattices are given in Table 1). 
The recent progress in experimental techniques em- 

ting optical tweezers llal. atomic force microscopy 
and soft microneedles [17| makes it possible to mon- 
itor the behavior of various polymers under tension and 
stress. In particular, applying a force on an isolated pro- 
tein, the mechanism of force-driven phase transitions was 
studied, such as unfolding of giant titine protein [lg| . 
the stretching and unzipping of collapsed DNA molecules 
fl9l ] . Of special interest in biophysics is the stretching of a 
collapsed polymer, i.e., of a polymer in a poor solvent be- 
low the 0-temperature. Unfolding proteins in this way 
could give important information on their spontaneous 
folding pathways. An intriguing question is whether in- 
termediate stages observed in different protein denatu- 
rating processes under applying force would always cor- 
respond to folding pathways in an enforced environment 
[20|.l2ll.l22l|. Note that varying temperature alone, a poly- 
mer cannot acquire the conformation of the completely 
stretched state with size exponent v = 1. Hence force 
not only influences the elastic, mechanical and struc- 
tural properties of polymers, but also introduces a new 
stretched state [regime of string) which is otherwise not 
accessible. The properties of force-induced globule-coil 



TABLE I: Values of connectivity constants and 0- 
temperatures for SASAWs on pure regular lattices and site- 
diluted percolative lattices for different space dimensions d. 
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and coil-string transitions in d = 2, 3 have been studied 
intensively both analytically l24l 12a. l26l H3, [H[ and ap- 
plying numerical simulations [26|, [29| • The existence of 
many intermediate states was found when unfolding un- 
der an applied force, which are not the same as the ther- 
modynamically stable intermediate states. The phase 
diagrams of the globule-coil and coil-string coexistence 
when varying temperature and stretching force were con- 
structed. 

New challenges have been raised recently in study- 
ing protein folding in the natural cellular environment 
[301 ] . Real biological cells can be described as a very 
crowded environment built of the biochemical species, 
which occupy a large fraction of the total volume; the 
"volume exclusion" arises due to the steric repulsion be- 
tween molecules [3l|. In the language of lattice models, 
the crowded environment with structural obstacles can 
be considered as a disordered lattice, where some amount 
of randomly chosen sites contains defects. Of particular 
interest is the case, when the concentration p of lattice 
sites allowed for the SAWs equals the critical concen- 
tration Pc (p c {d=2) = 0.592746 [32], p c (d=3) = 0.31160 
(HI) and the lattice becomes percolative. An incipient 
percolation cluster appears in the system, on which in- 
finitely long SAWs can reside. In this regime, the scaling 
law for the end-to-end distance holds with a new expo- 
nent v Vc > ^saw [HBSSSimiillii A hint to 

the physical understanding of this phenomenon is given 
by the fact that the percolation cluster itself is a frac- 
tal object with fractal dimension d^ a dependent on d. In 
this way, scaling law exponents of SAWs change with the 
dimension d^ of the (fractal) lattice on which the walk 
resides. Note that studying the scaling of SAWs on a per- 
colative lattice, one is interested rather in the backbone 
of the percolation cluster, which is defined as follows. 
Assume that each site of the cluster is a resistor, the 
neighbour sites are connected by conducting bonds and 
an external potential drop is applied at two ends of the 
cluster. The backbone is the subset of the cluster consist- 
ing of all sites through which the current flows; i.e., it is 
the structure left when all "dangling ends" are eliminated 
from the cluster. The SAWs can be trapped in "dangling 
ends" , therefore infinitely long chains can only exist on 
the backbone of the cluster, which is characterized by its 
own fractal dimension d£ (d^(d = 2) = 1.650 ± 0.005 
HH, d£(d = 3) = 1.86 ±0.01 0). 

Whereas the scaling behavior of SAWs on percolative 
lattices served as a subject of numerous studies since 
the early 80th (see, e.g., Ref. [13] for a recent review), 
less attention has been paid to peculiarities of the in- 
fluence of the fractal structure of the underlying lattice 
on properties of the coil-globule transition. The upper 
critical dimension shifts to d c = 6 for both SAW and 
0-point statistics on a disordered lattice at the percola- 
tion threshold. The scaling of polymer size at the coil- 
globule transition point is governed by > vq for 
d < 6 (e.g., vg{d = 2) = 0.74 ± 0.02 g|, vg(d = 
3) = 0.60 ± 0.02 (HI). It is established that the value 



of the 6-temperature is lowered due to the presence of 
disorder g3, |H, 0, EE H| ES HI , but estimates for 
were found up to now only for the case of bond percola- 
tion. The existing estimates are: TQ°(d = 2) = 0.62(6) 
0, Tg c (d = 3) = 0.43(6) [H (note, that corresponding 
values of connectivity constant for SAW on bond-diluted 
percolative lattices read: /j, Pc ~ 1.29 and 1.05 for d = 2 
and 3, respectively [42J). 

The response of a polymer in crowded media to the 
stretching force modelled by the SASAW model on per- 
colative lattice has been considered so far only in d = 2 
for relatively short chains by exact enumeration [49] ]. 
However, much more important is studying this problem 
in three dimensions, which describes real polymer sys- 
tems. This still needs a careful analysis and clarification. 
Also, a quantitative description of the globule-coil tran- 
sition under applying force, in particular the estimate of 
0-temperatures in disordered environment under stretch- 
ing, still remains an open question. 

The aim of the present study is to apply state-of- 
the-art numerical simulations to analyze the peculari- 
ties of SASAWs on site-diluted lattices at the percola- 
tion threshold (modelling flexible polymers in crowded 
environment) under applied external stretching force in 
space dimensions d = 2, 3. We estimate the shift of 
the 0-temperature of the globule-coil transition under 
the influence of stretching and analyze the effect of ap- 
plied force on the phase transitions between collapsed, 
extended and stretched phases. 

The outline of the rest of the paper is as follows. In the 
next section we describe the details of the numerical al- 
gorithm used in our study. In section III we present our 
results of coil-globule transition peculiarities for SAWs 
on percolative lattices, and in section IV we analyze the 
influence of stretching force on properties of the transi- 
tion between collapsed and extended states. We end up 
by giving conclusions and an outlook in section V. 



II. THE METHOD 

We consider site percolation on regular lattices of edge 
lengths up to L max =400, 200 in dimensions d=2,3, re- 
spectively. Each site of the lattice was assigned to be oc- 
cupied with probability p c and empty otherwise. To ob- 
tain the backbone of a percolation cluster on a given dis- 
ordered lattice, we apply an algorithm consisting of the 
following two steps: first finding the percolation cluster 
based on the site-labeling method of Hoshen and Kopel- 
man [50( and then extracting the backbone of this cluster 
[5l| (the al gori thm is explained in detail in our previous 
papers [iol |52| ) . We constructed 1000 clusters in each 
space dimension. 

To study SASAWs on the backbone of percolation clus- 
ters, we apply the pruned-enriched Rosenbluth method 
(PERM) [ll|, taking into account that a SASAW can 
have its steps only on the sites belonging to the back- 
bone of the percolation cluster. PERM is based on the 
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original Rosenbluth-Rosenbluth (RR) method [53| and 
enrichment strategies [54| . The polymer grows step by 
step, i.e., the nth monomer is placed at a randomly cho- 
sen empty neighbor site of the last placed (n — l)th 
monomer (n < N, where N is the total length of the 
chain). The growth is stopped, if the total length of the 
chain is reached. In order to obtain correct statistics, 
any attempt to place a monomer at an already occupied 
site would result in discarding the entire chain. This 
leads to an exponential "attrition" (the number of dis- 
carded chains grows exponentially with the chain length) . 
The bias due to avoiding this case is corrected in the RR 
algorithm by means of giving a weight W n ~ 11^=2 mi 
to each sample configuration at the nth step, where 
mi is the number of free lattice sites to place the Ith 
monomer. This method is particularly useful for studying 
0-polymers, since the Rosenbluth weights of the statisti- 
cally relevant chains approximately cancel against their 
Boltzmann probability. 

Population control in PERM suggests pruning config- 
urations with too small weights, and enriching the sam- 
ple with copies of high- weight configurations . These 
copies are made while the chain is growing, and continue 
to grow independently of each other. Pruning and en- 
richment are performed by choosing thresholds W„ and 
W„ depending on the estimate of the partition sums of 
the n-monomer chain. If the current weight W n of an 
n- monomer chain is less than W^, a random number 
r=0, 1 is chosen; if r=0, the chain is discarded, otherwise 
it is kept and its weight is doubled. Thus, low- weight 
chains are pruned with probability 1/2. If W n exceeds 
W„ , the configuration is doubled and the weight of each 
copy is taken as half the original weight. For updating 
the threshold values we apply similar rules as in JH^, H(| : 
W r >=C(Z„/Z 1 )(c„/c 1 ) 2 and W<=Q.2W>, where c„ de- 
notes the number of created chains having length n, and 
the parameter C controls the pruning-enrichment statis- 
tics. After a certain number of chains of total length 
N is produced, the iteration is finished and a new tour 
starts. We adjust the pruning-enrichment control param- 




FIG. 1: Self-avoding walks on the backbone of a percolation 
cluster in d — 2 under stretching force F. 



eter such that on average 10 chains of total length N are 
generated per each iteration 56], and perform 10 6 itera- 
tions. 

One end of the chain is subjected to an external force 
F acting in a chosen direction, say x (see Fig. [T]), while 
the other end (the starting point) is kept fixed. The 
stretching energy E s arising due to the applied force for 
an n-step trajectory is given by: 



E, 



(1) 



where x is the x-component of distance from the starting 
point \x n — xq\. The Rosenbluth weight factor W n is thus 
taken to be: 



n 

1=2 



-(S l -B,_ 1 )+F|x l -x,_ 1 | 



m ; e 



(2) 



where Ei denotes the energy of the /-step chain (Ei = Zi-e 
with e being an attractive energy between two nearest 
neighbors and Z\ the number of nearest neighbors con- 
tacts for a given chain) and ks is the Boltzmann con- 
stant. In what follows, we will assume units in which 
e = -1, k B = 1. 

When a chain of total length N is constructed, a new 
one starts from the same starting point, until the desired 
number of chain configurations is obtained. The config- 
urational averaging, e.g., for stretching in x-direction, is 
then given by 
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where Xg ont and 



'N 



lf denote ^-coordinates of start- and 
end-points, respectively, and Wff nl is the weight of an 
TV-monomer chain in a given configuration. 

Note, that studying SAWs on disordered lattices, we 
have to perform two types of averaging: the first average 
(...) is performed over all SAW configurations on a single 
backbone of percolation cluster; the second average (...) 
is carried out over different realizations of disorder, i.e., 
over many backbone configurations. 



III. 9-TRANSITION OF SASAW ON 
PERCOLATIVE LATTICE 

Statistical fluctuations of the energy E of a polymer 
chain, expressed by the behaviour of the specific heat 
Cy, signalize thermodynamic activity in the system, and 
thus the peak structure of Cy as a function of temper- 
ature indicates transitions or crossovers between physi- 
cally different states. In the case of a polymer system, 
this corresponds to the transition between globule and 
coil regimes. Cy can be expressed via energy fluctua- 
tions as follows: 



Cy(T) 



1 

2^2 



(E 2 ) - (E)- 



(4) 
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FIG. 2: Left: Energy distribution for N = 180-step SA- 
SAWs on a pure lattice in d = 2 at different temperatures: 
T = 0.6 (squares), T = 1.0 (diamonds), T = 1.5 (triangles). 
The broadening of the distribution curve emerges close to 
T e = 1.499 ±0.002. Right: Energy distribution for N = 180- 
step SASAWs on the backbone of percolation clusters in d — 2 
at different temperatures: T = 0.2 (squares), T = 0.4 (dia- 
monds), T = 0.9 (triangles). The broadening of the distribu- 
tion curve emerges close to TK a = 0.92 ± 0.02. 



It is worthwile first to discuss the behaviour of the en- 
ergy distribution Pr(E) with varying temperature of the 
system, as presented in Fig. [5] Since the energy dis- 
tribution shows one peak only, the transition could be 
denoted as being second-order like. The width of the dis- 
tribution grows with increasing temperature, until it has 
reached its maximum broadening value. This happens in 
the vicinity of the O-point. At higher temperatures, the 
distribution becomes narrower again. 

For comparison with available results and testing our 
methods, we performed also simulations on pure lattices. 
Figures [3j [4] show the specific- heat behavior of SASAWs 
on a pure lattice and the backbone of percolation clusters 
in d = 2, 3 for several different chain lengths. Note, that 
the maximum of the specific heat per monomer grows 
with N for SAWs on a pure lattice, whereas for the case 
of the backbone of percolation clusters it descreases with 
increasing N. For finite chain length TV, the temperature 
defined by position of specific heat maximum T™^ X (N) is 
well below the collapse transition O-temperature. This 
finite-size deviation of T™*(N) from Tq obeys scaling 



FIG. 4: Specific heat per monomer as a function of tem- 
perature for a SASAW in d = 3 on a pure lattice (left) and 
backbone of percolation cluster (right). Squares: N = 40, 
diamonds: N = 80, triangles: N = 120. 



behavior with N. Fig. [5]presents the chain-length depen- 
dence of the specific-heat peaks for pure and percolative 
lattices. For the case of pure lattice, d — 3 is the upper 
critical dimension for 0-transition, and the mean-field 
Flory-Huggins theory [57[ suggests: 



rpmax 
1 G V 



(5) 



In the general case, the approach of T™ X (N) to its lim- 
iting value obeys: 



T™y X (N) - T e ~ a ■ N~ v& + 



N 



(6) 



where a, b are constants, u& is the size exponent of a SAW 
at O-point: 

d = 2, pure lattice: v@ = 4/7 @, (7) 
d= 2, pc: i/ e c = 0.74 ±0.02 0, (8) 
d = 3, pc: i/ e c = 0.60 ±0.02 0. (9) 



Note, that we consider the special case of collapse tran- 
sition on the site-diluted percolative lattice, which was 
not studied before. 

Our estimates for Xq c are obtained by least-square fit- 
ting of d6j) , taking into account (HJ) and ([9]), and pre- 
sented in Table 1. As it was already mentioned, the 
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FIG. 3: Specific heat per monomer as a function of tempera- 
ture for a SAW in d — 2 on a pure lattice (left) and backbone 
of percolation cluster (right). Squares: N = 50, diamonds: 
N = 100, triangles: N = 150. 



FIG. 5: Peak temperatures of the specific heat of a SASAW 
in d = 2 (left) and d = 3 (right) at different chain length N. 
Squares: pure lattice, diamondes: backbone of percolation 
cluster. 
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FIG. 6: Left: Energy distribution for N = 180-step SASAWs 
on pure d = 2 lattice at T = 0.4 under a force F acting in the 
environment: F = 0.1 (squares), F = 1.0 (diamonds), F = 
1.5 (triangles). Right: Energy distribution for N = 180-step 
SASAWs on the backbone of percolation clusters at T = 0.1 
under force acting in the environment: F = 0.1 (squares), 
F = 1.0 (diamonds), F = 1.5 (triangles). 



collapse transition temperature is a non-universal prop- 
erty, which depends on the lattice type and, in partic- 
ular, on the connectivity constant. As expected, the 
values of Tq appear to be smaller than the pure lat- 
tice values, as far as relation fi Pc < /i holds. In addi- 
tion, due to the fact that /tt Pc decreases with d, the same 
tendency is reflected in the transition temperature be- 
haviour: 7* c (d=2) > Tg c (d=3). This is in contrast to 
the values on the pure lattice, but coincides with the ten- 
dency of bond-diluted percolation values, cited in the In- 
troduction. Let us note, finally, that our values for Tq of 
site percolation are larger than that of bond-percolation, 
which can be explained again by the difference of connec- 
tivity constants of these two fractal lattice structures. 



FIG. 8: Specific heat per monomer of a SASAW with N = 90 
steps in d — 3 under stretching force F as a function of tem- 
perature. Left: pure lattice, right: backbone of percolation 
cluster. Squares: F = 0.2, diamonds: F — 0.4, triangles: 
F = 0.6. 



(see Fig. 1). The energy distributions Pr(E) at fixed 
temperature and varying F are presented in Fig. [6] for 
the cases of the pure lattice and the backbone of perco- 
lation clusters for comparison. We have chosen the value 
of T in both cases to be well below the corresponding 
0-temperatures, so that at F = we restore the energy 
distribution of the globular state. With increasing F, 
the averaged energy of the chain decreases - applied force 
stretches the polymer globule. The width of the distribu- 
tion changes at increasing F, and reaches its maximum 
broadening in the vicinity of the O-point. The value of 
the transition temperature is now shifted by the presence 
of force. At higher values of F, the distribution becomes 
narrower again - the chain is an in extended state. 
To study the 8-transition of SASAWs, when the ex- 



IV. 9-TRANSITION OF SASAW ON 
PERCOLATIVE LATTICE UNDER 
STRETCHING FORCE 

Let now act an additional stretching force F in the re- 
direction of SASAW trajectory, fixed at its starting point 
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FIG. 7: Specific heat per monomer of a SASAW with N = 90 
steps in d — 2 under stretching force F as a function of tem- 
perature. Left: pure lattice, right: backbone of percolation 
cluster. Squares: F — 0.2, diamonds: F = 0.4, triangles: 
F = 0.6. 



FIG. 9: Peak temperatures of the specific heat of a SASAW 
as a function of chain length N. Upper row: d — 2, left: 
pure lattice, right: backbone of percolation cluster. Lower 
row: d = 3, left: pure lattice, right: backbone of percolation 
cluster. Squares: F = 0, diamonds: F — 0.2, triangles: F = 
0.4, pluses: F = 0.6. 
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FIG. 10: Phase diagrams of the stretching of a SASAW under 
applied force F. Upper row: d = 2, left: pure lattice, right: 
backbone of percolation cluster. The solid line presents re- 
sults of exact enumeration 29]. Lower row: d = 3, left: pure 
lattice, right: backbone of percolation cluster. The filled tri- 
angle present result of [1| F = 1.5, Te » 2.46. 




FIG. 11: Upper row: Averaged extentions of a SASAW under 
applying force acting in a d — 2 environment at T — 0.2; left 
- pure lattice, right - backbone of percolation cluster, squares: 
N = 20, diamonds: N = 50, triangles: N = 100. Lower row: 
Averaged extentions of a SASAW under applying force acting 
in a d — 3 environment; left - pure lattice (T — 1.8), right - 
backbone of percolation cluster (T = 0.4), squares: N = 20, 
diamonds: TV = 50, triangles: TV = 100. 



ternal stretching force is acting in the environment, we 
are working in the "constant-force" ensemble. Fixing the 
value of F, we study the specific-heat behavior (Figs. [3 
[8]). Analyzing the peak structure of the specific heat, 
we immediately conclude, that increasing the value of F 
leads to decreasing the transition temperature. Fig. [5] 
presents the chain-length dependence of the specific-heat 
peaks for a pure and percolative lattices at several differ- 
ent values of F. 

For finite chain length N, the temperature defined by 
the position of the specific-heat maximum T™ X (N) is 
well below the collapse transition O-temperature. Our 
estimates for Tq in the presence of force are obtained 
by least-square fitting of © with 0, ©, @. Results 
are presented in Fig. [TU] in the form of a phase diagram 
of transitions from globule to extended state, and listed 
in Table 2. For the pure lattice in d = 2, we compare 
our results with an exact enumeration study of Kumar 
et al. [29|, where SASAWs of length up to N = 55 under 
stretching force were studied. Our results also appear to 
be in good correspondence with that of Ref . [23] ■ 



TABLE II: Values of O-temperature on the backbone of per- 
colation clusters in d — 2, 3 at varying force F. 



F 




0.2 0.4 0.6 0.8 1.0 


rpPc 

1 e 

rpp c 
2 & 


(d = 
(d = 


2) 0.92(2) 0.84(2) 0.74(3) 0.56(2) 0.40(4) 0.18(5) 

3) 0.71(2) 0.66(2) 0.60(2) 0.48(3) 0.35(2) 0.20(2) 



Averaged extention in ^-direction, giving information 
about the internal structure of the polymer configuration 
under applying force, is presented in Fig. [Tl]for the cases 
of pure lattice and backbone of percolation cluster. At 
small forces, a polymer chain is still in the compact folded 
state and slightly oriented along the force direction. At 
larger forces, the polymer chain has a conformation sim- 
ilar to the extended (swollen) structure. Note, that com- 
pletely stretched states, corresponding to (x)/N ~ 1 can 
be obtained only on the pure case and are not accessi- 
ble on the percolative lattices due to complicated frac- 
tal structure on the underlying percolation cluster. Our 
estimates of (x) /N in d = 2 appear to be in a good 
correspondence with that of available exact enumeration 
studies HI,!!!]. 



V. CONCLUSIONS 

We studied the self-attracting self-avoding walks on 
disordered lattices in space dimensions d = 2, 3, mod- 
elling flexible polymer macromolecules in porous environ- 
ment. We considered the special case, when the concen- 
tration of disorder is exactly at the percolation threshold, 
so that an incipient percolation cluster of sites, allowed 
for SAWs, emerges on the lattice. In our study, SASAWs 
reside only on the backbone of the percolation cluster, 
which has a fractal structure. 

In the first part of our study, attention has been paid 
to the influence of the fractal structure of the underlying 
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lattice on properties of the coil-globule transition. Ap- 
plying the pruned-enrichcd Roscnbluth algorithm, we ob- 
tain estimates of the collapse transition temperature Xq c 
of SASAWs on site-diluted percolative lattices in d = 2 
and 3 dimensions (note, that so far only estimates for 
the bond-percolation case have been found). The values 
of Xq c appear to be smaller than the pure lattice value. 
This can be explained, remembering that the collapse 
transition temperature is proportional to the connectiv- 
ity constant fi of SAWs on a given lattice. In addition, 
due to the fact that [i Pa decreases with d, the same ten- 
dency is reflected in the anomalous transition tempera- 
ture behaviour: Tg c (d=2) > Tg c (d=3). 

Next, keeping one end of a SASAW trajectory on 
the backbone of percolation cluster fixed, we applied a 
stretching force, acting in some chosen direction (say, 
x). Especially interesting was to study this problem in 



d = 3, which corresponds to real polymer systems and 
was not considered before. We estimated the shift of the 
0-temperature of the globule-coil transition under the 
stretching and constructed phase diagrams of collapsed 
and extended states coexistence in d — 2, 3. As expected, 
the presence of stretching force in environment leads to 
a decreasing 0-temperature value. 
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